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Abstract 

The stability of He@C6o and Ne@C6o is discussed in the context of a spherical 
model where the carbon atoms are smeared out into a uniform shell. The 
electronic properties of the sixty 7r electrons together with those of the central 
atom are treated in the Thomas-Fermi approximation. Simple electrostatic 
reasoning elucidates the nature of the radial stability of the complex. A 
method to include non-spherical corrections is outlined. Possible bonding 
topologies of the central atom and the Ceo cage are discussed, as well as the 
relevance of these topologies to incipient central atom distortions. 
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I. INTRODUCTION 



An interesting consequence of the closed cage structure of C 60 is that atoms and molecules 
may be trapped inside, forming endohedral complexes ( "dopeyballs" ). A variety of such com- 
plexes have been produced [0 where the dopant may be a metallic or a rare gas atom/cluster. 
The electronic structure and mechanical stability of two such endohedral rare gas complexes, 
He@C6o and Ne@C6o, are considered here. 

It has previously been noted P] that the high symmetry of the Ceo cage suggests the 
following geometrical approximation: the one-electron potential of icosahedral symmetry 
is replaced with its spherical average. Physically, the nuclear charge together with the 
charge from the core electrons is smeared out into a thin spherical shell of uniform surface 
charge density and radius R. Such a fictitious molecule- "spherene" -has complete spherical 
symmetry. 

In the case of metallofullerenes such as Na@C6o, it has been shown || that the equilibrium 
position of the dopant is not at the center of the cage, consequently reducing the symmetry 
of the complex. For rare gas dopants, however, it has previously been asserted [|l| that there 
is no such symmetry reduction. Additionally, previous calculations [[| on Ne@C6o indicate 
that the dopant is stable at the center. Consequently, given the large coordination number of 
the gas atom and the lack of a distortion from the high symmetry configuration, endohedral 
rare gas fullerenes seem to be good candidates for the spherical approximation. 

The endohedrally-doped spherene is treated in the Thomas-Fermi approximation. It is 
anticipated that the highly delocalized 7r-electrons of the cage and the closed shell config- 
uration of the gas atom are well-suited to such a statistical description. This method of 
treating high symmetry molecules was originally suggested by March ||. 

The cage is parametrized by the radius R and valence electron number N, while the gas 
atom has nuclear charge Z and an equal number of electrons. The total electron density is 
calculated using the Thomas-Fermi equations, subject to the shell boundary conditions. 

The total energy of the system is calculated for a variety of radii. A stable equilibrium 
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radius is obtained, when the cage self-energy is computed from the system of ionic point 
charges, rather than from the continuum model. A method of treating non-spherical cor- 
rections is outlined. Lastly, possible central atom distortions are discussed in the context of 
the Bader molecular graphs arising from the electron density. 



II. THOMAS-FERMI THEORY OF DOPEYBALLS 

The nuclear potential from the central atom and the cage atoms can be expanded in a 
multipole series as 



Qi m is the 2 £ -pole moment, given by 
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where Z is the atomic number of the rare gas atom of interest. From the icosahedral 
symmetry of the molecule, it has been shown || that the first three non- vanishing multipole 
moments are for £ =0, 6, 10. As the occupied one-electron states have an effective £ = 5, it 
is not necessary to consider £ > 10 multipole terms in first order perturbation theory. 

The spherical approximation consists of retaining only the £ = term in the multi- 
pole expansion. The error introduced can be estimated by consideration of the relevant 
dimensionless parameters 
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The following values are found as function of Z: a 6i0 = 0.026/ (1 + Z), a 6t5 = 0.020/ (1 + Z), 
a w,o — 0.021/ (1+Z), and a Wi5 = 0.034/(l + Z). As ai m <C 1 for £ < 10, it may be concluded 
that the spherical approximation is justified for such a high symmetry structure, and one- 
electron splittings under the true icosahedral symmetry can be treated perturbationally. 

With spherical symmetry imposed, the endohedral cluster is quasi-atomic in form. The 
sixty 7r-electrons of the cage plus the Z electrons of the central atom are now treated in the 



Thomas- Fermi (TF) approximation where the Z + N electrons see a point charge of Ze at 
the origin and a uniform shell of charge Ne and radius R. The methodology follows that of 
March [||. 

At temperature T = 0, the dimensionless TF equation without exchange is given by 
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x is the distance from the center of the shell in units of 
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where ao is the Bohr radius of hydrogen. is related to the potential in the usual way 

V(r) = —Mx) 
r 

Eq. [|is supplemented with the boundary conditions: (f)(0) = 1, and 

N 
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where X is the shell radius in dimensionless units and differentiation is with respect to x. 
Additionally, itself is continuous over its domain, and — > as x — > oo. 



III. NUMERICAL RESULTS 

0(r) is found by numerical integration of Eq. ^ subject to the above boundary conditions. 
The charge density n(r) and total electronic energy E e are subsequently obtained from 0(r), 



as 
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as was obtained by March. The radial electron density of Ne@Ceo is given as an example in 
Fig. § 

The electronic energy E e is found for various values of the dimensionless cage radius 
X in the cases of Z = 2 and 10, corresponding to He@Ceo and Ne@Ceo respectively. As 
the boundary condition for at the shell explicitly depends on Z, E e will not have simple 
"atomic" scaling with Z (E e ~ Z 7 ^ 3 ). 

If we add to E e the electrostatic self-energy of the continuum shell, U, it is found that 
this total energy E = E e + U does not have a minimum for finite X, in accord with Teller's 
no-binding theorem || for molecules in TF theory. The continuum shell self-energy is too 
large and dominates E e . However, if the self-energy is evaluated as a sum over point ions 
distributed on the shell surface, an energy minimum at finite shell radius is found. The 
spherical approximation is abandoned for the purposes of calculating the potential energy 
for the nuclear configuration. 

The self-energy of the nuclear configuration may be written in the form 

„ ZNe' iVV 

c/ = ^r +c i3r (10) 

where c is a dimensionless number, computed from the actual equilibrium coordinates of 
Ceo- It was previously |2j computed as c = 0.43101. 

The total energies for He@C 60 and Ne@C 60 as a function of the shell radius are displayed 
in Figs. |] and [| The resulting equilibrium radii, R , are listed in Table |, together with 
the total energies for the equilibrium configuration. Previous results for Cgo are provided 
for comparison. 

The expansion of the equilibrium cage radius with increasing Z is observed, in agreement 
with a previous restricted Hartree-Fock calculation [Q. It should, however, be noted that 
c is a result of the detailed atomic positions in Cqq. Certainly, distortions which do not 
preserve the relative positions of cage ions will give rise to changes in c, and consequently, 
to changes in Rq. It is also interesting to note that the monotonic increase of Rq with Z 
does not follow the "atomic size" scaling relation (D ~ Z^ 1 ^ 3 ). 



The shell equilibrium can be understood from electrostatic considerations. At equilib- 
rium, the centrifugal force on the shell resulting from self-interaction is balanced by the 
centripetal force exerted on the shell by the charge contained inside the shell. The total 
charge contained inside the shell at equilibrium is —cNe ~ — 25.86e, and the total number 
of electrons inside the shell is (Z + cN). Thus we see that endohedral fullerenes whose 
central atoms have higher Z require more electrons inside the shell to maintain equilibrium. 



While approximating the discrete cage ions by a uniform spherical shell is valid approx- 
imation with regard to the total energy, it is the corrections to the spherical average which 
contain all the information regarding the bonding. The spherical approximation reduces the 
complex to a one-dimensional system. The consequent critical points in the electron den- 
sity consist of only minima and maxima, and an analysis of the bonding topology requires 
consideration of the effects of the nuclear multipole moments £ > 0. A method for finding 
the electron density corrugations of endohedral fullerenes is given below. 

To treat these non-spherical contributions, one must return to the general TF equation, 



IV. BEYOND SPHERENE 
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where f3 = 32 3 ^ 3 6 (2m) 3//2 and p + is the (positive) nuclear charge density. 



p + can be expressed as a sum of spherical harmonics, 
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where 




(13) 



The potential V is now written as 
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where Vo is the spherically-averaged solution and ( results from consideration of the higher 
order multipole moments. It is assumed that ( is much smaller than V . Thus, Eq. [XT] leads 
to a linearized equation for ( 

V 2 C = - AnAp (15) 

Eq. is of the form of the single particle Schrodinger Eq. at zero energy with a central 
"potential," ^f3V , which is everywhere positive, and an non-homogeneous boundary term. 
C may be expanded in spherical harmonics inside and outside r = R, and the boundary term 
gives rise to a discontinuity in the radial derivative of (. Thus, 

((rVEWW) (16) 
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where Q m satisfies the following 
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and is subject to the following boundary conditions 



r&m->0, as r^O (18) 
rQm^O, as r oo (19) 
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As the effective "potential" in Eq. [T7| is positive for all r, only exponential solutions 
are possible for £. Without a discontinuity in the radial derivative, only the trivial solution 
would satisfy the boundary conditions at the origin and at infinity. Thus, (e m is non-zero only 
when Qi m is non-zero. The first non-spherical corrections to the potential and the electron 
density then are at t = 6, 10. Not surprisingly, these £ values correspond to irreducible 
representations of the rotation group which contain the trivial representation (A g ) of Ih- It 
is interesting to note that the centrifugal term in the effective "potential" will reduce the 
contributions from large I to the density. 
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V. STRUCTURE AND BONDING 



The simplest bonding configuration consistent with the symmetry constraints is that of 
sixty bonds between the central atom and the carbon atoms in the cage. However, a second 
possibility was proposed on the basis of an ab-initio calculation |D on Ne@Ceo- There, it 
was found that with the Ne atom in the center of the cage, thirty bond paths exist, starting 
on the Ne atom and terminating on the carbon-carbon double bond points. Is this a stable 
configuration with respect to central atom displacement? 

While the total energy calculations in Ref. H indicate that Ne is in a stable equilibrium 
at the center of the cage, unless the bond point at the carbon-carbon double bond is a non- 
nuclear attractor |§ ( "pseudoatom" ) , it would seem that the bonding topology found in 
Ref. || may correspond to that of the metastable state described by Bader as a "conflict 
structure." 

Given that the Na + is isoelectronic with Ne, and that Na + in Na + @Ceo does distort 
along a five-fold axis, a similiar instability in Ne@Ceo would seem possible. The possible 
incompatibility of the molecular graph with the total energy calculations indicates that 
additional study is warranted, and central atom distortions in rare gas endohedral fullerenes 
remain as an intriguing possibility. 
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FIGURES 




FIG. 1. Radial electron density for Ne@Ceo (in units of 10/6) vs x for R = 7.36ao- 
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FIG. 2. Total energy E (Ry) for He@C 60 vs shell radius R (Bohr). 
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FIG. 3. Total energy E (Ry) for Ne@C 60 vs shell radius R (Bohr). 
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TABLES 



TABLE I. Minimum total energies E (Ry), equilibrium radii R (Bohr), for Z = 0, 2, 10. 
Z RE 
a 7.35 -25.313 

2 7.41 -32.98 

10 7.64 -356.30 

a 55 . i = 
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